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We propose orbital fluctuations in a multi-band ground state as the superconducting pairing 
mechanism in the new iron-based materials. We develop a general SU(4) theoretical framework for 
studying a two-orbital model and discuss a number of scenarios that may be operational within 
this orbital fluctuation paradigm. The orbital and spin symmetry of the superconducting order 
parameter is argued to be highly non-universal and dependent on the details of the underlying band 
structure. We introduce a minimal two-orbital model for the Fe-pnictides characterized by non- 
degenerate orbitals that strongly mix with each other. They correspond to the iron u d xy " orbital 
and to an effective combination of "dzx" and ll d zy " , respectively. Using this effective model we 
perform RPA calculations of susceptibilities and effective pairing interactions. We find that spin 
and orbital fluctuations are, generally, strongly coupled and we identify the parameters that control 
this coupling as well as the relative strength of various channels. 

PACS numbers: 



I. INTRODUCTION 

The recently discoveredi^i2i^'ii&£ iron-based super- 
conducting oxides with a transition temperature as high 
as 55K bring up the immediate important question of 
the mechanism underlying superconductivity in this new 
class of materials. An obvious temptation is to connect 
the superconductivity in these new materials to that in 
the high-T c cuprates based on a number of compelling 
similarities: layered 2D nature of both classes of com- 
pounds, multi-element oxide nature of the materials in 
both cases, key role of doping with the undoped mate- 
rials being non-superconducting, importance of nearby 
magnetic states whose suppression leads to superconduc- 
tivity, relatively high superconducting transition temper- 
atures. We argue here theoretically, using very general 
considerations, that in spite of these tantalizing similar- 
ities between the cuprates and the new Fe-based super- 
conductors, there is a very important qualitative differ- 
ence which suggests that the nature of superconductiv- 
ity in the Fe-based oxides is likely to be different from 
that of the cuprates. In particular, the hallmark of the 
Fe-based superconductors seems to be the multiband na- 
ture of their low energy band structure as revealed by 
first principles band structure calculation !] 9 ' 10 ! 11 ! 12 ! 13 ] 14 , 
and as such, orbital fluctuations are likely to be a crucial 
ingredient of physics in these new materials. We propose 
in this article that the superconductivity in the Fe-based 
oxides is driven by novel orbital fluctuations (which are 
invariably coupled to the spin fluctuations, as described 
below) in the multi-band ground state which have no ana- 
log in the standard model of the high-T c cuprates. The 
existence of multiple bands at the Fermi surface gives rise 
to a new paradigm, namely, the possibility of a supercon- 
ducting pairing mechanism driven by orbital fluctuations, 
which we believe are the underlying cause for the high- 
temperature superconductivity in these materials. 

To explore the possibilities opened by this new 
paradigm we study an effective two-band tight-binding 



model of the FeAs planes. Band structure calcula- 
tions show that the electronic character near the Fermi 
level is mostly determined by the d-orbitals of Fe~ and 
that all five orbitals participate in the formation of the 
Fermi surface^. However, there are several reasons that 
make the study of a simple two-band effective model 
highly relevant: Firstly, a two-band model is the ideal 
minimal model that combines the requirements for ob- 
serving multi-orbital physics with a relative simplicity 
that makes this physics transparent. Secondly, LDA- 
type calculation o 9 ' 10 ! 11 ' 12 ! 13 ' 14 have shown that the un- 
dope d 16 i 17 ' 18 i 19 material is characterized by five Fermi 
sheets: two quasi-2D electron pockets located near the 
M point of the Brillouin zone and two quasi-2D and one 
3D hole pockets in the vicinity of the T point. How- 
ever, upon electron doping the 3D hole pocket disap- 
pears and the 2D hole pockets shrink rapidly. More- 
over, strong-coupling LDA+DMFT calculations^ show 
that at finite doping the bands responsible for the elec- 
tron pockets clearly cross the Fermi level, while the hole 
pockets around the F point are barely identifiable. These 
preliminary results suggest that the low-energy physics 
responsible for superconductivity in the iron-based com- 
pounds may be governed by the two-bands associated 
with the electron pockets and, therefore, an effective two- 
band model would be an appropriate way to describe 
it. Nonetheless, as we show below, the nature of the or- 
bitals that participate in the formation of these bands, as 
well as the strength and symmetry of the hybridization 
between them have direct consequences for the pairing 
mechanism. Furthermore, to completely clarify the key 
question of the relative importance of the electron and 
hole pockets in the low-energy physics of the FeAs lay- 
ers as a function of doping, further calculations involving 
momentum-dependent self-energy effects are necessary. 

The general analysis described in the next section has 
three main objectives. First, in order to disentangle and 
describe unambiguously various possible types of fluctua- 
tions we introduce a set of generators of the SU(4) group 
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associated with on-site rotations. We identify the oper- 
ators that describe charge, spin as well as charge-orbital 
and spin-orbital coupled degrees of freedom. We then ex- 
press the bare interaction, as well as the effective interac- 
tion in terms of these operators and identify the couplings 
corresponding to the relevant channels. The second ob- 
jective is to classify all possible pair operators, which in 
this particular case are combinations of spin and orbital 
singlets and triplets, and express the interaction in terms 
of these pair operators in a manner suitable for mean- field 
calculations. The coupling constants corresponding to 
each possible pairing channel will be expressed in terms 
of the effective interactions introduced in step one. We 
identify the channels susceptible to generate pairing and 
determine the role of orbital fluctuations in the mech- 
anism leading to superconductivity. We show that the 
bare interactions can not generate pairing, but the inclu- 
sion of orbital and/or spin fluctuations opens multiple 
possibilities for pairing. The specific spin-orbital sym- 
metry of the superconducting order parameter is highly 
non-universal and depends on the details of the underly- 
ing band structure and on the coupling constants. 

The paper is organized as follows: In Section [IT] we de- 
scribe in detail a general framework suitable to describe 
the physics of the SU(4) spin-orbital space that charac- 
terizes a two-orbital effective model. We also introduce a 
distinction between inter-orbital and inter-band pairing 
and discuss several general pairing scenarios that may 
be at work in a two-orbital system, emphasising the role 
of the orbital degrees of freedom. In Section IIIII we in- 
troduce our effective two-band model for the Fe-pnictides 
and compare it with other two-orbital models. Using this 
effective model, we calculate the spin and orbital suscep- 
tibilities and the corresponding renormalized interactions 
at the random phase approximation (RPA) level. The 
results are presented in Section [iVl We show explicitly 
the role of the orbital degrees of freedom and point out 
the highly non-universal nature of the pairing problem in 
multi-orbital systems. Section IVl contains a summary of 
this work and our final conclusions. 



II. GENERAL SU(4) FORMALISM FOR 
TWO-ORBITAL MODELS 



We start with a two-orbital minimal model given by 
the Hamiltonian H = H t + Hjj. Explicitly, the kinetic 
term is 
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are hopping matrix elements between the n 



where VJf 1 

orbital at site j and the m orbital at site i, tr 6 {1\ 1} IS 
the spin label and /i is the chemical potential. The en- 
ergy bands for the non-interacting system are obtained by 
diagonalizing the matrix e„ m (k), the Fourier transform 
of tYi n . The interaction term contains contributions from 
the on-site intra-band (U) and inter-band (If) repulsion, 



as well as the Hund's coupling J and pair hopping J', 

Hu = U E n im] n iml + ^' E n ^ n ^' 
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Next, we introduce the generators X^ v of the SU(4) 
group defined as2£ 



— ^ ] ^ ] c ims °ss' r mm' c im's't (3) 
s.s' m.m' 

where a = ^2x2 is the identity matrix and a a with 
a 6 {1, 2, 3} are Pauli matrices and for a two-band model 
and f M = <t m . There is a total of 16 X-operators. No- 
tice that X®° represents the charge, while = 2Sf 
represents the pure spin. Similarly, T" = X® a are the 
charge-orbital operators, while the nine operators X"^ 
with a, (3 € {1,2,3} correspond to the spin-orbital cou- 
pled degrees of freedom. 

Including the effects of fluctuations, the effective inter- 
action Hamiltonian has the form 



Hf" = ^ 
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where Xq are Fourier components of the X-operator and 
the factor 1/16 was introduced for later convenience. In 
the most general case T is a 16 x 16 matrix. However, 
as we are interested in studying fluctuations of a param- 
agnetic state characterized by SU(2) spin symmetry, the 
effective interaction Hamiltonian is rotationally invariant 
in spin space and the matrix T becomes block-diagonal, 
with one 4x4 block corresponding to the charge and 
charge-orbital degrees of freedom and three 4x4 iden- 
tical blocks corresponding to the spin and spin-orbital 
coupled degrees of freedom. Explicitly we have 

fpM p'M) = r$(q) S p0 5 p/0 + rW(q) 5 pa 5 p , a , (5) 

where a £ {1,2,3}. The effective coupling matrices 
r( c )(q) and r( s )(q) can be obtained diagrammatically 
starting with the bare couplings. Before writing explic- 
itly these bare couplings in terms of the original inter- 
action parameters U, U', J and J', let us make two 
observations. First, we note that the two-particle op- 
erators appearing in the interaction Hamiltonian ^ are 
generated by diagonal combinations X^X^, while the 
off-diagonal contributions X^X^ with fj, ^ v generate 
new operators that do not appear in Hjj, so we expect 
the matrices r°( c ' s ) to be diagonal. Second, we notice 
that, while expressing Hjj in terms of X-operators is al- 
ways possible (up to a constant energy shift), there is 
no unique choice for the coupling matrices r°( c / s '. For 
example, terms proportional to the on-site intra-orbital 
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interaction operator D = J2 m c lmT c ™J- are generated in 
the charge channel, Xf®X®° , the charge-orbital channel 
X? 3 X? 3 , as well as in the spin channel X?°X? , and the 
spin-orbital channel X" 3 X" 3 . Consequently, there are 
distinct linear combinations of products of X-operators 
equal to D. Our choice is based on the following two 
conditions: i) the linear combinations of products of X- 
operators should be invariant under spin rotations, and 
ii) to each channel that generate a certain term of the 
interacting Hamiltonian we will ascribe a proportional 
fraction of the corresponding coupling constant. Using 
this procedure we obtain 

r°( c ) = diag{(£7 + 2U' - J), (-[/' + 2 J + J'), 
{—U' + 2J— J'), (U-2U' + J)} 

T ^ = dmg{(-U-J), (-U'-J% (6) 
{-U' + J'), (-U + J)} 

We begin our second task by classifying the pair oper- 
ators A a f,(k) defined as linear combinations of products 
c Ls c -km's" We choose combinations that are symmet- 
ric or antisymmetric under permutations of the spin (or- 
bital) indices, representing the spin (orbital) triplet and 
singlet channels, respectively. In the notation A ab (k) 
the first index refers to the orbital degree of freedom 
and takes the value a = s for the orbital singlet and 
a £ {t\,t2,to} for the orbital triplets. The spin index is 
b = s for the spin singlet and b £ {£+, <o} for the spin 
triplets. With these notation the 16 independent pair 
operators are 
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A a b (k) = - [s r cT 1T cL k2i + s 2 ci cli c T _ k2T 

+ S 3 c lk2T C -kl| + 54 C k21 C -klT 



where m £ {1,2}, a £ {T,i} and the set of signs 
(si, S2, S3, S4) takes the values (+,+,+,+) for (ab) — 
(t t ), (+,-,+,-) for (ab) = (t s), (+,+,-,-) for 
(ab) = (sto), and (+,—,—,+) for (ab) = (ss). Note that 
all the singlet-singlet and triplet-triplet pair operators 
are odd functions of momentum, A ss (— k) = — A ss (k) 
and A tt (— k) = — A tt (k), consequently pairing involving 
these operators occurs in the p-wave channel. On the 
other hand, the singlet-triplet and triplet-singlet pair op- 
erators are even functions of k, A st (— k) = A st (k) and 
A ts (— k) = A is (k), and are involved in s-wave or d-wave 
pairing. 

Using Eq. ([7]) we express the effective Hamiltonian ^ 

as 

H u f *T, E W aba ,»(ky)A ab (k)Al b ,(k'), (8) 

k,k' a,b : a' : b' 



where, as usual, we retain only contributions of the 
form cl _ „ cJ „ , „ cu' m ' s > m > s > with k 2 = — ki and 
k' 2 = — k' x . Because of the symmetry properties of 

the pair operators, most of the elements of the W ma- 
trix are in fact zero. The non- vanishing elements are 
W ss — W ss ss describing the pairing interaction in the 



orbital-singlet spin-singlet channel, W s t = W t 



ss t b t b : 



with 



tb £ {£+,£_, to}) representing the orbital-singlet spin- 
triplet pairing, plus two 3x3 matrices [Wtt]t t a = 

w t ai n t a2 t„ and [W ts }t ai t a2 = W tai s t a2 s corresponding 
to the orbital-triplet spin-triplet and orbital-triplet spin- 
singlet channels, respectively. Note the the coupling ma- 
trices are degenerate with respect to the spin-triplet label 

Next we express the pairing interaction W ab in terms 



(c/s) 



In order to 



of the effective coupling constants T 
simplify the expressions, we will take into consideration 
only the diagonal contributions T^ c / S ^ and T^ s ^ = 
which is always non-zero for non-degenerate orbitals, and 
assume that all the other coupling constants are negligi- 
ble. With these assumptions, we have for the singlet- 
singlet and singlet-triplet channels 



Wst = ^ [U 
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triplet-triplet pairing matrix becomes 
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A similar expression can be written for the triplet-singlet 
channel, Wt s , by making the substitution T+, — > T~. 
Expressing the pairing interaction in terms of the effec- 
tive coupling constants T ^ completes our second objec- 
tive and provides the tools necessary for understanding 
orbital fluctuations mediated pairing. 

Using the bare couplings T° given by Eq. ©, the 
pairing interaction (|9l fT0)) becomes 



000 
000 
000 



wP = 4 



W» = 0; 







U' + J 



W° t = -(U'-J) 



(11) 



At this level, pairing could only appear in the s-wave 
channel and have cither orbital-triplet spin-singlet or 
orbital-singlet spin-triplet character. The necessary con- 
ditions are U' < J, or U < J' . If we use parameters 
characterizing isolated Fe atoms, these conditions are cer- 
tainly not satisfied, as J, J' <C U,U'. Even if we consider 
that the two-band effective model should contain renor- 
malized values of these couplings, it is not very likely 
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that these conditions, especially U < J', be realized™ 
We conclude that, at the bare level, pairing is unlikely 
to occur. Theoretically, it may be realized only in the 
orbital-singlet spin-triplet channel if U' < J 22 . The sta- 
bility of this channel in the presence of fluctuations will 
be discussed in section HVl 

Before discussing the role of fluctuations, we need 
to clarify the concepts of inter-orbital and inter-band 
pairing. In the absence of hybridization between or- 
bitals, tfj n — (or e nm (k) = 0) for n =^ m, the wave- 
vectors k and — k characterizing pair operators of the 
type Cjj ncr c_ kmo ./ cannot be chosen to lie on the Fermi sur- 
faces corresponding to the different orbitals, except when 
they are degenerate or at points of accidental degen- 
eracies 2 ^^. The inter-pocket orbital-singlet spin-triplet 
pairing discussed in Ref. corresponds to this situa- 
tion. For the sake of clarity, we will use for this case the 
term inter-band pairing. In contrast, for a non-zero hy- 
bridization, the relevant momenta k and — k should be 
considered on the Fermi surface produced by the energy 
bands E n (k) obtained by diagonalizing Ht (not on the or- 
bital "Fermi surfaces" generated by e nn (k)). Assuming 
inversion symmetry, the pair created by such an inter- 
orbital pair operator is well defined for any Fermi momen- 
tum. Note that diagonalizing the hopping Hamiltonian 
H t is equivalent to performing a unitary transformation 
u(k). Using this transformation we can express orbital 
creation operators in terms of band creation operators 
as c\. ma = ttm^(k)<4 M(r . Consequently, we can project 
the orbital pair operator into the low energy band and 
obtain a band pair operator relevant for the low-energy 
physics. Note that the two operators may have differ- 
ent symmetries^ 2 -. To set the stage for our discussion of 
the orbital fluctuations, we conclude that, because the 
phase space available for inter-band pairing is in general 
very limited, this type of pairing is not likely to play a 
significant role in the mechanism responsible for super- 
conductivity in the iron-based compounds. In contrast, 
inter- orbital pairing has no phase-space limitations and 
could be one of the key elements that characterizes the 
physics of these multi-orbital superconductors. 

In the presence of hopping, the effective interac- 
ts^ 

tions r^V will acquire a momentum dependence. For 
example, at the RPA level, the bare couplings © 
will be renormalized by the susceptibility X/^(q) 
■ Em,)!!',!),!! 1 C„ m (k + q)T mm ,G m i n i (k)r^, n , where 
G„ m (k) is the Green function ((ckn CT cj cmCT )). It is conve- 
nient to express the elements of the susceptibility matrix 
in terms of Xmm', nn'(q) =EkC™™'( k + q)G™'(k). The 
diagonal elements become 

X00(q) = ^ (XllAl + Xl2,21 + X21.12 + X22,22) , 
Xll( c l) — 2 (^11,22 + X12.12 + X21.21 + X22,ll) , (12) 
X22(q) = ~ (Xll,22 — Xl2,12 — X21.21 + X22,ll) , 



X33(q) = ^ (Xll,H — Xl 2 , 2 l — X 2 1,1 2 + X22,22) ■ 

The most significant off-diagonal components are 
Xo3(30) = (Xii,ii±Xi2,2iTX2i,i2 -X22,22)/2, which 
couple the spin susceptibility xoo and the longitudinal 
orbital-spin component X33 ■ Note that the rest of the off- 
diagonal components x^v depend on contributions of the 
form xa,ij or Xijjj that are typically much smaller than 
the diagonal terms. Neglecting, for simplicity, the off- 
diagonal susceptibilities, the re-normalized coupling con- 

stants become r^ /s) (q) = T°^ /s) / (l + r° ( ; /s) Xwi (q)) . 

Based on this expression, together with equations ([6]), 
(fT2|) and ([9l fT0|) we identify three simple pairing scenar- 
ios. To establish the exact relation between these scenar- 
ios and the physics of the iron-based superconductors, 
detailed numerical calculations are necessary. Also, for 
certain values of the band and interaction parameters, 
more complicated situations that do not involve a cer- 
tain dominant component of the susceptibility are possi- 
ble (see Section HV|) . 

Scenario A. Xnn,nn ^ > {Xnn,mmi Xnm,n'm'} i i.e., the 

intra-orbital susceptibilities are dominant. This can be 
realized, for example, when the Fermi surface is charac- 
terized by an approximate nesting with wave-vector Q 
and the hybridization is negligible. The dominant ele- 
ments of the susceptibility matrix are xoo and X33 ■ Con- 
sidering the expression © of the bare interaction, wc 
conclude that the strongest renormalization will occur in 

(s) 

the spin channel, T 0Q , followed by the spin-orbital chan- 

(s) . (c) 

nel Tgg and the charge-orbital channel T 33 . However, in 
the presence of any non-zero Hund's coupling J, the pure 
spin channel will dominate. This translates into a large 
negative contribution to Fq (Q) or a large positive renor- 
malization of r^~ (Q). Regardless of channel, spin fluctu- 
ations represent the driving force in this scenario. How- 
ever, because X03 is generally non-zero, spin fluctuations 
couple to the longitudinal spin-orbital fluctuations, thus 

(s) 

enhancing T 33 . Very importantly, this coupling between 
the spin and spin-orbital channels generates stronger fluc- 
tuations and, consequently, a stronger renormalization 
of the effective coupling constants in both channels. No- 
tice that orbital fluctuations, which renormalize r^ 3 , give 
contributions to W st and W ss , as well as to W to t and 
Wt o8 , which are opposite in sign to the contributions from 
the spin channel and, consequently, are detrimental to 
pairing. In contrast, for Wt a t and Wt aS , spin and orbital 
fluctuations have a similar effect on pairing. Consider- 
ing also that in T~ the term r( s ) has a coefficient —3, 
while in T + the coefficient is +1, we conclude that the 
most likely pairing in this scenario is orbital-triplet spin- 
singlet pairing. This implies s-wave or d-wave symmetry 
for the order parameter. 

Scenario B. {Xnn,nniXnm,7nn} Xnn,mmi i.e., the 

intra-orbital components are large and the hybridization 
is strong. In this case the off-diagonal terms x nm ,mn 
contributing to xoo and X33 become important. The sign 
of these contributions at the relevant wave-vector Q de- 



pends on the symmetry of the hybridization ei2(k). A 
positive contribution will further enhance the spin fluc- 
tuations. On the contrary a negative contribution makes 
X33 larger that xoo and this could offset the effect of the 
Hund's coupling J and promote the spin-orbital coupling 

(s) 

T 33 as the main component of the pairing interaction. 
Pairing in this scenario is mainly mediated by either or- 
bital fluctuations or spin fluctuations, depending on the 
sign of Xnm.mm i- e -j on the symmetry of the hybridiza- 
tion. Note that the key ingredient leading to the further 
enhancement of spin and spin-orbital fluctuations is the 
strong orbital hybridization. We emphasize that, in ad- 
dition to the strength, a key role is played by the sym- 
metry of the inter-orbital hybridization, as it determines 
the sign of the off-diagonal susceptibilities Xnm,mn- As in 
scenario A, the coupling constants that are most strongly 

(s) (s) 

renormalized are T QQ and Igg , leading to orbital-triplet 
spin-singlet pairing as the most likely type of paring, with 
s- or d-wave symmetry for the order parameter. 

Scenario C Xnn,mm \Xnn,nn: Xnm.n'm'} : i.e., the 

inter-orbital components are larger than the intra-orbital 
contributions in the vicinity of certain relevant Q vec- 
tors. In this scenario the elements xu an d X22 of the 
susceptibility matrix become the dominant components. 
Consequently, the main contribution to the renormal- 
ized pairing interaction will come from the effective cou- 
plings Tn and ]?22- The pairing mechanism is now con- 
trolled by orbital fluctuations, as in scenario B for neg- 
ative off-diagonal susceptibilities Xnm,mn- However, this 
time it is the inter-orbital susceptibility that leads to a 
strong increase of the effective coupling in the charge- 
orbital and especially the spin-orbital channels. This 
will translate into a positive momentum-dependent con- 
tribution to the even-parity channels W st and W toS or a 
negative momentum-dependent contribution to the odd- 
parity channels W ss and W to t- However, because in T~ 

the rjuw contribution comes with a pre- factor -3, we can 
infer that always singlet-singlet or the triplet to - singlet 
channels will be preferred over the triplet to -triplet and 
the singlet-triplet channels, respectively. 



III. EFFECTIVE TWO-ORBITAL MODEL 

For a specific discussion of the orbital effects in Fe- 
pnictides, we need a simple tight-binding model that 
describes these materials. In particular, a choice of 
the hopping parameters in the non-interacting Hamil- 
tonian (fTJ) has to be made. In making this choice, 
we take into account the results of band structure 
calculation a 9 ' 10 ! 11 ' 12 ! 13 ' 14 and try to reproduce as well as 
possible the low-energy sector. The goal is to reproduce 
not only the multi-pocket structure of the Fermi sur- 
face, but also the correct energy scales for the low-energy 
modes. We note here that effective two-orbital models 
have been already discussed in the literatur o 25 i 26 ' 27 i 28 ' 29 , 
but they focus exclusively on the degenerate d zx , d zy or- 
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FIG. 1: (color online) : Energy dispersion for the effective two 
orbital model along the (0, 0) — > (tt, 0) — * (tt, tt) — > (0, 0) path 
in momentum space. Orbital "1" corresponds to the iron 
d xy orbital, while "2" corresponds to an effective d zx - d zy 
combination. The dispersion of the second orbital is shown 
for a set of hopping parameters corresponding to the presence 
of relatively large hole pocket (see main text), i5if 2 )(k) is 
the diagonalized high (low) energy band. Inset: Zero energy 
contours inside the unfolded Brillouin zone. The F point of 
the original Brillouin zone corresponds to {[0, 0], [tt, tt]}, while 
M corresponds to {[tt, 0], [0,7r]}. 

bitals. However, d xy is also known to play an impor- 
tant role in the formation of the electron pockets and it 
strongly mixes with the two degenerate orbitals. Poten- 
tially important aspects of pairing mechanism may occur 
because of this mixing between the d xy orbital and cither 
d zx or d zy , suggesting a three-orbital minimal model for 
the Fe-pnictides^ 2 -. Our goal is to capture some of this 
physics within a two-orbital model. Consequently, we 
consider an effective model of the Fe-based oxides consist- 
ing of two non-degenerate orbitals: the first corresponds 
to the iron d xy orbital, while the second represents an "ef- 
fective orbital" that simulates the combination of d zx and 
d zy . More precisely, the energy band associated with the 
second orbital approximates the low energy band struc- 
ture of d zx hybridized with d zy . 

The intra-orbital hopping parameters that define 
the model are t%™ = (t[o,o] , <[i,o] , £[1,1] , £[0,2]), where 
[Ax, Ay] = [{x,j — Xi),{yj — yi)] denotes a given hop- 
ping vector, as well as those related to it by symmetry. 
The first orbital is characterized by the set of hoppings 
t\j = (0.2,-0.32,0.17,-0.015) eV. The proper energy 
scale is set by the nearest-neighbor hopping in i = —0.32 
e V 13 i 30 . For the second orbital we use three different sets 
of parameters. If tff = (t 22 ) Q = (0.18,0.22,0.16,0.07) 
eV, a hole pocket is generated around the T = (tt, tt) 
point (see Fig. []}. Again, the energy scale for the dis- 
persion curve is determined by the choice of the nearest- 
neighbor hopping in accordance with the band structure 
calculation a 13 ! 30 . For a smaller hole pocket, which would 
correspond to the doped case, we use t 22 = (tfj)„ = 

(0.18, 0.22, 0.14, 0.05) eV. To simulate the case when the 
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FIG. 2: (color online): Comparison of the band structure of 
our effective two-orbital model {E\, E2) with the band struc- 
ture of the d zx - d zy model from Refs. HHH]] (E+, E-) 
corresponding to a nearest-neighbor hopping ^1=1 eV and 
a chemical potential fi — 1.45 eV. Notice the similarities be- 
tween the high-energy bands E\ and E+ and the different 
energy scales of the lower energy bands. 

hole pockets are absent, or have a very weak coherent 
spectral weight^, we use the set of hopping parameters 
t?| = = (0.18,0.22,0.125,0.04) eV. Note that the 

only significant change in the dispersion curve occurs in 
the vicinity of the T point, where the hole pocket is lo- 
cated (for details, see the inset of Fig. [3}. Finally, the 
hybridization between the two orbitals is described by 
£12 (k) = 4(5 sin fc^ s\ak y with 5 = 0.25eV. The dispersion 
curves, diagonalized bands and zero energy contours are 
shown in Fig. [TJ Note that our model reproduces the 
two electron Fermi pockets predicted by LDA calcula- 
tions (in Fig. [TJ M corresponds to the [0, tt], [n, 0] points 
of the unfolded Brillouin zone), but generates only one 
hole pocket at V (i.e., near [tt, it}). We expect the absence 
of a second hole pocket to play a minor role in the qual- 
itative aspects of the pairing mechanism, as one pocket 
is enough for capturing the physics stemming from the 
approximate nesting between the electron and the hole 
Fermi pockets. 

To place our effective model in the context of similar 
two-band models for the iron-based superconductors, we 
show in Fi 
from Refs 



ig. El a comparison with the d z 
. I2I2I2I The high energy bands, E 1 (k) and 



d Z y model 



E + (k), respectively, which are responsible for the forma- 
tion of the electron pockets have the same energy scale 
and very similar dispersions. By contrast, the low-energy 
band E_(k) of the d zx — d zy model has a characteristic 
energy that is an order of magnitude larger than the low- 
energy band E^fo) of the present model. Particularly 
significant is the rapid dispersion of E-(k) in the vicin- 
ity of [7r,7r], which results, on the one hand, in a much 
larger Fermi energy for the hole pocket and, on the other 
hand, in smaller values of the susceptibility. Our model 
correctly takes into account the results of band struc- 
ture calculations that give characteristic Fermi energies 
of about 0.2 eV for both types of pockets. 




(0,0) (7C,0) (7C,Jt) (0,0) (71,0) (71,71) (0,0) 



FIG. 3: (color online): Left panel: static spin susceptibility 
(Xoo-red line) and static transverse orbital susceptibility (xn- 
green line) for the two-orbital model studied in Ref. [26|. The 
dash line represents the maximal eigenvalue of the matrix. 
Right panel: largest components of the spin-orbital suscepti- 
bility matrix in the RPA approximation. Notice that both 
Xoo an d X11 diverge near (ir/2,iv/2), signaling the presence 
of an instability. The instability is produced by the coupling 
of the spin (pi = 0) and transverse spin-orbital (fj, = 1) modes, 
due to a non- vanishing xoi • 

IV. RPA SUSCEPTIBILITIES AND EFFECTIVE 
PAIRING INTERACTIONS 



The random phase approximation (RPA) involves 
the summation of the whole series of bubble diagrams 
and, within the formalism described in Section [Til 
it corresponds to the multiplication with the matrix 

(l + f °( c ' a 'x(q)) • The result depends critically on 
the structure of the matrices f°( c / s ) and %(q). Let 
us consider the case corresponding to our model (and 
also to the degenerate orbital model), when the inter- 
orbital hybridization is real and has even parity, i.e., 
ei 2 (— k) = ei2(k). Within the SU(4) formalism, the cou- 
pling constant matrix is diagonal, with components given 
by Eq. ©, while the bare susceptibility matrix has the 
form 



(13) 



/ Xoo 


Xoi 





X03 


X10 


Xu 





Xia 








X22 





V X30 


X31 





X33 



where the nonvanishing components are 
1 



Xoo (33) 

Xll (22) 
X03 

X01 

X31 



:(Xll,ll + X22,22 i 2xi2,2l); 



Xll, 22 ± Xl2,21) 

Xll, 12 
Xll, 12 



X30 

Xw 

Xl3 



X22,22), 



X22,21, 
X22,21- 



(14) 



(15) 



Note that the transverse orbital component X22 is decou- 
pled from the other channels and its RPA value will de- 
pend only on T®^^ ■ By contrast, the spin (charge) sus- 
ceptibility xoo always couples to the longitudinal orbital 
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0.2 1 1 1 . 1 , 1 

(0,0) (Ji,0) (7i,7i) (0,0) 

FIG. 4: (color online): Intra-orbital susceptibility ^22,22 corre- 
sponding to three different sets of hopping parameters (fff) , 
x £ {a,/3,7}, for the second orbital (see main text). The cor- 
responding dispersion curves in the vicinity of F (i.e., [it, it]) 
are shown in the inset and are characterized by (a) a large 
hole pocket (red squares), (/3) small hole pocket (green cir- 
cles), or (7) no hole pocket (blue diamonds). The existence 
of the hole pocket determines the appearance of peaks in the 
susceptibility near [tt,0], due to an approximate nesting be- 
tween the electron and the hole Fermi surfaces. 

component X33- I n addition, in the presence of an inter- 
orbital hybridization, these two components are coupled 
to the transverse orbital susceptibility xn- All these cou- 
plings can significantly enhance the susceptibilities. To 
give an example, we show in Fig. [3]the static spin suscep- 
tibility Xoo^ — 0) and the transverse orbital susceptibil- 
ity Xn( w — 0) f° r the two-band model studied in Ref. 
|26| . The maximal eigenvalue of the susceptibility matrix 
(dashed line in Fig. [3]) is peaked near q = [ir/2,n/2] 
as a results of the mixing between xoo and Xn- Note 
that Xoi(io)(q) ^ while X3i(i3)(q) = X30(03)(q) = 
because of symmetry. At the RPA level, for an interac- 
tion U = U' = 2.8 and J = J' = both xoo and xn 
are strongly peaked near q (right panel in Fig. [3]), in- 
dicating the proximity of a spin-orbital instability. This 
instability is the result of the the spin channel (// = 0) 
being coupled to one of the transverse spin-orbital chan- 
nels (/j, = 1) due to the hybridization between orbitals. 
The strength of this coupling is determined by the off- 
diagonal susceptibility Xoi(io)- Note that the numerical 
value of Xoi(io) is negligible everywhere except in the 
vicinity of q where it has a peak Xoi(io) ~ —0.1. 

Next, we turn our attention to the RPA analysis of 
the effective two-band model described in the previous 
section. Using this simple model we want to address two 
basic questions: i) what is the dominant pairing channel, 
and ii) what is the specific contribution of the orbital 
degrees of freedom. We find that the physics of the two- 
band model is rather non-universal, depending strongly 
on both the structure of the non-interacting Hamiltonian 
and the form of the interaction. However, within a rea- 
sonable parameter window we find that (i) the dominant 
pairing channel is the orbital-triplet spin-singlet chan- 
nel, and that (ii) the spin and orbital components are 




(0,0) (71,0) (jt.TC) (0,0) (71,0) (Jt,7t) (0,0) 



FIG. 5: (color online): RPA susceptibilities in the spin sec- 
tor for various parameters of the two-orbital model: a) Band 
structure with a small hole FS and no Hund's coupling, b) 
Model with a small hole pocket and large Hund's coupling, 
c) No hole pocket and J — 0. d) No hole pocket and large 
Hund's coupling. Note that the mixing between the pure spin 
and the longitudinal orbital-spin components, measured by 
X03, is particularly strong in the absence of a strong Hund's 
coupling. 

strongly coupled and typically give comparable contribu- 
tions to the effective pairing interaction W. The rela- 
tive weight of the spin channel is enhanced by increas- 
ing the size of the hole pocket, which leads to nesting, 
and by increasing the Hund's coupling J. In the opposite 
limit the orbital effects become dominant. In addition, 
we note that the symmetry of the hybridization, which 
determines the off-diagonal susceptibility X1221, plays a 
crucial role in determining the relative strength of various 
contributions. Our present choice, ei2(k + Q) = — £12 (k) 
for Q = (it, 0) or Q = (0,7r), enhances the spin compo- 
nent over the orbital-spin contributions. By contrast, an 
even symmetry choice would further enhance the role of 
orbital fluctuations. We show that these coupled spin- 
orbital fluctuations are important and cannot be ignored 
in intrinsically multi-orbital systems like the iron-based 
superconducting oxides. In terms of the general scenar- 
ios described at the end of section [TTl we find scenario 
B as the most likely to be realized within our effective 
two-orbital model. 

The first step in the RPA analysis involves the calcu- 
lation of the bare susceptibilities XijA'j' f° r each of the 
three sets of parameters describing our the two-orbital 
model, which were discussed in section IIII1 The most 
significant parameter dependence is seen in X22.22 as the 
second orbital is responsible for the existence and size 
of the hole pocket. We show in Fig. [4] the intra-orbital 
susceptibility X22,22 corresponding to the three different 
sets of hopping parameters (if?) , x € {a,(3,j}. The 
case (a) corresponds to a large hole pocket (see the inset 
of Fig. [4]) that produces an approximate nesting with 
the electron pocket and generates a strong peak at [it, 0]. 
This features becomes weaker as the hole pocket shrinks 
- case (/?) - and eventually disappears in the absence of 
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the hole pocket - case (7) - being replaced by a wide 
maximum near [tt, tt]. 

What happens in the presence of an interaction? To 
reduce the number of independent parameters we only 
consider the case 

U' = U- 2 J, J' = J. 

In the absence of Hund's coupling, J = 0, the RPA 
spin channel susceptibility x^Hq) diverges at the crit- 
ical values of the interaction Uc"^ ~ 0.75, Uc® « 0.85 
and Uc^ ~ 1.2 for the three sets of orbital parameters, 
respectively. In the presence of a hole pocket the insta- 
bility occurs in the vicinity of [tt, 0], while if only the 
electron pockets exist, the instability occurs at [7r,7r]. In 
both cases there is a strong mixing between the pure 
spin and the longitudinal orbital-spin components while 
the transverse orbital-spin components are practically de- 
coupled. Including a Hund's coupling J strengthens the 
[tt, 0] instability or replaces the orbital-driven [tt, tt] in- 
stability with a spin-driven [tt, 0] instability. The general 
trends of the RPA susceptibility are shown in Fig. [5] for 
band structures corresponding to two sets of band pa- 
rameters, (3 (small hole pocket) and 7 (no hole pocket), 
and interactions characterized by either J = or by 
a relatively strong Hund's coupling, J/U ~ 0.35. We 
note that a non-vanishing J always enhances the spin 
susceptibility, while suppressing the longitudinal orbital- 
spin component. However, in a more realistic model ex- 
change interactions between nearest-neighbor and next- 
nearest-neighbor sites have to be considere d 17 ' 18 ' 19 ' 21 . In 
this case, the bare coupling constant matrices r°( c / s ) will 
acquire momentum dependent contributions J(q). De- 
pending on the character (ferromagnetic, antiferromag- 
netic or mixed) of the exchange couplings, the spin sus- 
ceptibility can be further enhanced or, by contrary, sup- 
pressed in favor of the orbital components. We stress 
that the inclusion of the exchange interaction is crucial, 
as it can qualitatively change the susceptibilities and, im- 
plicitly, the nature of the possible instabilities and that 
of the pairing mechanism. 

We turn now our attention to the calculation of the 
effective pairing interaction W a b given by Eqns. © and 
(flTjf . We show in Fig. [6] the momentum dependence 
of the dominant effective couplings. The largest compo- 
nent is always the second diagonal element of the orbital- 
triplet spin-singlet matrix, [Wt s ] 22 . hi the absence of a 
hole pocket, for J = (panel (c) in Fig. [6]), [W / ts ] 22 (k) 
has a wide maximum at (tt, tt), while the other effective 
couplings have a much weaker momentum dependence. 
By contrast, if the hole pocket exists (panels (a) and 
(b)), [Wt s ] 2 2 (k) is characterized by sharp maxima in the 
vicinity of (tt, 0). These maxima are the consequence 
of the approximate nesting between the hole and elec- 
tron Fermi surfaces. Note that the orbital-triplet spin- 
singlet pair operator ||7J) is an even function of momen- 
tum. Consequently, because of the of the structure near 
(tt, 0), the superconducting gap will have s-type symme- 
try with opposite signs on the hole and electron pockets. 




(0,0) (71,0) (71, Jt) (0,0) (71,0) (71, je) (0,0) 



FIG. 6: (color online): Momentum dependence of the dom- 
inant components of the effective coupling for the parameter 
regimes of Fig. [5] Note that the largest effective interaction 
occurs in the orbital-triplet spin-singlet channel (red line) and, 
in the presence of a hole pocket - (a) and (b) - it leads to s- 
wave pairing with opposite signs of the gap on the electron 
and hole Fermi surfaces. The orbital dominated case (c) is 
consistent with d-wave pairing, while in case (d) there are 
several competing channels. 

By contrast, in the absence of the hole pockets and for 
J = 0, the (tt,tt) maximum will enforce a d-type sym- 
metry. The experimental evidence for the gap symme- 
try is rather contradictory. Some point-contact Andreev 
reflection experiment o 31 ' 32 , as well as microwave^ 3 - and 
ARPES 34 ' 35 measurements suggest fully-gaped supercon- 
ductivity, which is consistent with s-wave pairing. On the 
other hand, the results of other point-contact Andreev re- 
flection experiments 3 ^, together with NMR results^ 7 - and 
H c i magnetization measurements 3 ^, suggest a nodal gap 
function, consistent with d- or p-wave pairing. Clearly, 
further work is necessary in order to clarify this problem. 
Note that, regardless of symmetry, if paring occurs in 
the orbital-triplet spin-singlet channel it originates in the 
d Z x — d zy "effective" orbital. However, because the off- 
diagonal elements [W( S ] 12 = [Wt s ] 21 are non-zero, pairing 
is also induced in the d xy orbital. Therefore, the result- 
ing pairing gap has two components which are character- 
ized by two different energy scales that can be obtained 
by numerically solving a matrix gap equation. The gap 
equation can be derived using standard procedures from 
a BCS-type mean field approximation of the effective 
Hamiltonian ([8]). Notice that for large Hund's couplings 
(panels (b) and (d) in Fig. [6]) the relation J > U 1 is 
satisfied and the coupling constant in the orbital-singlet 
spin-triplet channel W s t becomes negative 22 , opening the 
possibility of s-wave pairing in this channel. Nonetheless, 
in the presence of a hole pocket - case (b) - because of 
the dominant contribution from [VF ts ] 22 , pairing realizes 
in the orbital-triplet spin-singlet channel. On the other 
hand, in case (d) - with no hole pocket and strong J 
- there are several closely competing channels. Finally, 
we note that for most of the parameter space our model 
predicts intra-orbital, rather than inter-orbital pairing. If 
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the Hund's coupling is small, scenario B is realized and 
pairing occurs in the orbital-triplet spin-singlet channel 
with A t2S as the strongest component. For large J, in 
the presence of a hole pocket, the relative importance of 
the off-diagonal susceptibility Xi2,2i diminishes and sce- 
nario A is realized. However, pairing still occurs in the 
orbital-triplet spin-singlet channel and is dominated by 
the second intra-orbital component. The only possibility 
for inter-orbital pairing occurs in the absence of a hole 
pocket and for large Hund's couplings, when pairing may 
occur in the orbital-singlet spin-singlet or in the orbital- 
singlet spin-triplet channels, depending of the details of 
the model. However, we emphasize again that these re- 
sults are highly non-universal and they strongly depend, 
for example, on the number of orbitals considered, on 
the symmetry and strength of the hybridization, or on 
the values of possible non-local interactions. 

V. SUMMARY AND CONCLUSIONS 

In summary, we propose a paradigm for superconduc- 
tivity in the new Fe-based materials where the pairing 
is caused by orbital fluctuations strongly coupled to spin 
fluctuations in the renormalized multiband ground state. 
We test our theoretical idea by carrying out an RPA-type 
calculation using a minimal two-orbital model consistent 
with the low energy bands determined by first princi- 
ples band structure calculations. We find that the intra- 



orbital pairing, rather than the inter-orbital or the inter- 
band pairings, mediated by coupled spin-orbital fluctu- 
ations is the driving mechanism here. One specific fal- 
sifiablc prediction of our paradigm is that, if the parent 
compound is not multiband (thus, rendering orbital fluc- 
tuations relatively unimportant), then the system would 
not be superconducting (or will have a rather low tran- 
sition temperature). Our numerical analysis of the the 
effective two-orbital model suggests several tasks and di- 
rections for future studies: i) It is crucial determine the 
size of the Fermi surfaces and the relative strength of 
the quasiparticles on the electron and hole pockets for a 
relevant doping range. If, for example, the quasiparti- 
cle residue takes significantly different values on the two 
types of Fermi surfaces, all the features coming from their 
approximate nesting will be strongly suppressed, ii) It is 
necessary to know exactly the symmetry of the inter- 
orbital hybridization and to have a good estimate of its 
strength. A model that just reproduces the correct low- 
energy band structure is not necessarily correct. The 
symmetry of the inter-orbital mixing terms plays a key 
role in establishing the relative weight and the couplings 
between various spin, charge and orbital modes, iii) A re- 
alistic tight-binding model of the Fe-pnictides should in- 
clude short-range interaction terms. The nature and the 
strength of these interactions have direct consequences 
on the possible instabilities, as well as on the supercon- 
ducting mechanism. 

This work was supported by LPS-NSA-CMTC. 
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